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The spin-boson model has nontrivial quantum phase transitions in the sub-Ohmic regime. For the 
bath spectra exponent ^ s < 1/2, the bosonic numerical renormalization group (BNRG) study 
of the exponents /? and 5 are hampered by the boson-state truncation, which leads to artificial 
interacting exponents instead of the correct Gaussian ones. In this paper, guided by a mean-field 
calculation, we study the order-parameter function m(r = a — Ofc, e. A) using BNRG. Scaling analysis 
with respect to the boson-state truncation Nb, the logarithmic discretization parameter A, and the 
tunneling strength A are carried out. Truncation-induced multiple-power behaviors are observed 
close to the critical point, with artificial values of P and 5. They cross over to classical behaviors 
with exponents /3 = 1/2 and 5 = 3 on the intermediate scales of r and e, respectively. We also find 
r/A^~^ and e/A scalings in the function m(r, e. A). The role of boson-state truncation as a scaling 
variable in the BNRG result for ^ s < 1/2 is identified and its interplay with the logarithmic 
discretization revealed. Relevance to the validity of quantum-to-classical mapping in other impurity 
models is discussed. 

PACS numbers: 05.10.Cc, 05.30.Jp, 64.70.Tg, 75.20.Hr 



I. INTRODUCTION 

The spin-boson model is the simplest model that de- 
scribes a quantum two-level system subjected to the in- 
fluence of a dissipative environment. It has applications 
in various fields in physics^ and its properties are stud- 
ied extensively;^ Especially, for the bath spectra exponent 
^ s < 1 (sub-Ohmic bath) , the ground state of the spin- 
boson model may change from the spin-tunneling state 
to the spin-pinned state through a second-order phase 
transition, as the dissipation strength crosses a critical 
value from below. This environment-induced quantum 
phase transition attracts much attention in the past few 
years.-"- 

Among the many theoretical methods that have been 
used to study this quantum phase transition, the bosonic 
numerical renormalization group (NRG) method is re- 
garded as the most accurate one. Technically, NRG is 
composed of three standard procedures: logarithmic dis- 
cretization, transforming the Hamiltonian into a semi- 
infinite chain, and the iterative diagonalization. Thanks 
to the logarithmic discretization, the state information at 
exponentially small energy scales is kept along the itera- 
tive diagonalization. Therefore, NRG allows for reliable 
extraction of critical exponents and description of the 
crossover behavior. Two parameters control the precision 
of NRG, i.e., the logarithmic discretization A (A > 1) 
and the number of kept states Ms. The original Wil- 
son's NRG is designed for the impurity problem with a 
fermionic bath.* In the past few years, the extension of 
NRG to impurity models with bosonic bath, such as the 
spin-boson model, proved to be fruitfuli^i^iiiS Due to the 
infinite local Hilbert space for each bosonic bath mode, 
one has to truncate the space into Ni, states. Usually 
the boson occupation states are used as the local bases, 
although optimal bases have also been consideredj^ Us- 



ing this bosonic NRG (BNRG), the thermodynamical as 
well as dynamical quantities for the spin-boson model 
are studied and the critical exponents in the sub-Ohmic 
regime obtained. Although it was clear that the localized 
fixed point (((Tz) ^ 0) cannot be described exactly due to 
the truncation of boson states, it was believed and partly 
checked^ that the BNRG shares the virtue of NRG, i.e, 
the parameters A, Mg, and Nt, only influences the value of 
non-universal quantities, such as the critical point value 
ttc and the prefactors of power laws. The exact value of 
them can be reliably obtained by extrapolating A ^ 1, 
Ms — > oo, and Nb oo. The universal quantities such 
as the critical exponents are not supposed to depend on 
A, Ms, or Nb. 

For the quantum critical behavior of the spin-boson 
model in the sub-Ohmic regime, the following critical ex- 
ponents have been studied^i^: 

m{a > (Xc, T = 0, e = 0) cx (a — 

m{a — ac, T — 0,e) (x e^^^ 

x(q!, T = 0, e = 0) oc \a — ad"'*' 

T*(a,e = 0) (X \a - 



x{a = ac,T,e ^ 0) (X T 



C{uj){a = ac, T 0, e = 0) cx cj-^ 



(1) 



Here, a and e describe the coupling strength between 
the spin and bosons, and the bias field on the spin, re- 
spectively. X C{uj) are static susceptibility and dy- 
namical spin correlation function, respectively. The naive 
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BNRG study of these critical exponents leads to the fol- 
lowing conclusions^jS (i) In the regime ^ s < 1, the 
critical fixed points are interacting and the correspond- 
ing critical exponents are non-classical; and (ii) The hy- 
perscaling relation and lo/T scaling hold. In the regime 
^ s < 1/2, these conclusions are in contrast to pre- 
vious theories based on the quantum-to-classical map- 
ping. There, the spin-boson model is mapped into a one- 
dimensional Ising model with J^- cx l/{ri — rj)(^+*).^ 
In the regime ^ s < 1/2, this Ising model is above 
its upper critical dimension, leading to Gaussian criti- 
cal fixed point and classical exponents. The hyperscal- 
ing relation and w/T scaling do not hold there j^i"— Re- 
cently, using a number of new methods,—^— the quantum 
phase transition in the sub-Ohmic spin-boson model has 
been studied. The obtained critical coup ling strength 
ac(s),— the exponent j/,— and 7 (Ref. are con- 
sistent with the BNRG results. However, the quantum 
Monte Carlo (QMC) simulationM and the exact diago- 
nalization study^- found that in the regime ^ s < 1/2, 
the critical point is Gaussian with classical exponent 
/3 — 1/2, being different from the BNRG conclusion. In 
the regime 1/2 ^ s < 1, BNRG results are consistent 
with the quantum-to-classical mapping theory which pre- 
dicts an interacting fixed point and oj/T scaling. 

Recently, a closer examination of BNRG method dis- 
closes two sources of error, which were not noticed before. 
One is the boson state truncation error^iS the other is 
the mass-flow error. The boson state truncation er- 
ror spoils the evaluation of the order-parameter related 
exponents /3 and S, while the mass-flow problem influ- 
ences the correct evaluation of x. The two sources of er- 
ror are different in nature and exist simultaneously in the 
BNRG algorithm in the whole regime ^ s < 1. But, 
they influence the critical behavior only in the regime 
^ s < 1/2, where the critical fixed point is expected to 
be Gaussian in the absence of these errors. 

For the mean-field spin-boson model which has a Gaus- 
sian critical fixed point, •^^ we showed that the boson state 
truncation leads to an artificial interacting fixed point 
in the regime ^ s < 1/2, but has no influence in 
1/2 ^ s < 1. This is an example where the boson-state 
truncation destroys the Gaussian fixed point and spoils 
the correct calculation of the exponents /3 and S. It leads 
to the surmise that the same may happen in the BNRG 
study of the spin-boson model. It would be difficult to 
find the Gaussian nature of the critical fixed point in 
^ s < 1/2 using BNRG, if the truncation works the 
same way as in the mean-field Hamiltonian. 

In Ref. |T^, the boson-state truncation error is traced 
back to the presence of a dangerously irrelevant variable 
for a Gaussian critical fixed point. The correct exponent 
can be seen on the intermediate scales. For the more fun- 
damental problem of mass- flow error, Vojta et al^ have 
proposed an extended NRG algorithm to partly remedy 
the problem and got the correct exponent a; = 1/2 in 
^ s < 1/2. In this paper, we focus on the boson state 
truncation error. In the BNRG, it is still unclear how a 



finite iVh leads to wrong /3 and S, and how to extract the 
correct exponents. For the spin-boson model, a thorough 
numerical study in the regime ^ s < 1/2 is required to 
prove or disprove the validity of the quantum-to-classical 
mapping in this model.— '^^i-'^^'^^i^^'^'^ Here, we use the 
scaling approach to analyze the BNRG data with respect 
to boson-state truncation Nh, logarithmic discretization 
parameter A, and tunneling strength A. We find that 
for any finite N^, the order parameter m has a multiple 
power form like that at the tricritical point ,"21 with non- 
classical exponents /3 and 6 dominated by the discretiza- 
tion scheme. The correct power- law behavior can be 
observed on the intermediate scale away from the crit- 
ical point. These two different power-law regimes are 
connected at a crossover scale, which goes to zero as a 
power of a; = 1/Ni, and w = A — 1. Thus, in the limit of 
either x ^ or w 0, the classical critical exponents /3 
and S are recovered. This is the same as in the mean-field 
spin-boson model, which we will detail in the Appendix. 
Besides, we also disclose the role of A as a scaling variable 
in the order parameter close to criticality. 

In Sec. II, the spin-boson model and our main results 
are presented. In Sec. Ill, a summary and discussion 
will be made. In the Appendix, we present the critical 
behavior of order parameter and susceptibility for the 
mean-field spin-boson model. 

II. MODEL AND RESULTS 

The Hamiltonian of the spin-boson model reads as 
H = -yCTj; + |(Tz + ^ (^ai + a}^ + ^uj,ala,. 

i i 

(2) 

Here, ax and CTz are Pauli matrices, and a.i and a] are the 
bosonic annihilation and creation operators of the mode 
i, respectively. The properties of the quantum two level 
system are determined by the environment through the 
following bath spectrum^: 

J{u;) = T:J2^1Siuj-u;,), (3) 

i 

for which we assume the following power form, 

J{uj) = 2TTauj'ujl-', (0 < a; < Wc). (4) 

Here, a measures the strength of the dissipation, cjc = 1 
is used as the energy unit. 

A. Mean-field results 

Before we carry out NRG calculations, it would be 
helpful to first have a look at the mean-field spin-boson 
model. It has Gaussian critical fixed point and classical 
exponents for any s ^ 0. It is used to mimic the situation 
of the full spin-boson model in the regime ^ s < 1/2. In 
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Ref. I21I the scaling behavior of to = (uz) with respect to 
boson state truncation Nj, was investigated numericahy. 
It was found that any finite N}, will lead to non-mean-field 
exponents /? and S in the regime < s < 1/2 and an ex- 
ponential behavior at s = 0. The exponents (as functions 
of s) agree well with those extracted from BNRG calcu- 
lations for the full spin-boson model. Here, we present a 
concise and complete analytical solution which will guide 
our BNRG study in the next section. 

The Hamiltonian of the mean-field spin-boson model 
reads as (neglecting a constant), 



H 



A 



i 

^(cr2)^Ai (a, -haj) +^a;,a|a,. (5) 



To make connection with the NRG study, we carry out 
a standard logarithmic discretization^iS. The obtained 
star-type mean-field Hamiltonian reads as 



TTStar 



2 ■ 



^ n 

^naian + '^0= ^ 7n + • (6) 



Here, the logarithmic discretization gives 
o 27ra 



7„ 



1 - A-(i+^) A-"(i+^)c^2^ 



and 



1 + s 1 - A-(2+«) 

2 + s 1 - A-(i+'5) 



A- 



(7) 



(8) 



A > 1 is the logarithmic discretization parameter. 

For the spin and boson decoupled Hamiltonian Eq.(6), 
the self-consistent equations are easily solved when there 
is no truncation, i.e., A^b = oo. For finite A^f,, this model 
cannot be solved exactly. However, through an analysis 
of the single-mode Hamiltonian, we get the asymptoti- 
cally exact expression for the order parameter to as a 
function of a — ac and e for given N^,^ A, and A. We 
summarize the results in the following and leave the de- 
tailed derivation in the Appendix. 

For a fixed A, the critical etc for Eq.(6) does not depend 
on A^b, but only on A. We define the following scaling 
variables r = a — ac(A), x — l/A'f,, and w = A — 1. To 
obtain the exponent we take e = 0. Magnetization 
m{T,e = 0,x,w) = {o'z)/2 has the following behavior in 
the limit t,x,w 0. For s ^ 1/2, 



For < s < 1/2, we get 
m(T, e = 0, a;, w) 

(f )^[a,(l)]^ (A) (^x)-^r^, (r « r,,); 



(r > Tcr). 



(10) 



Here, ac(l) = As/(2a;c) is the mean-field critical point 
for A = 1. The crossover scale Tcr reads as 

2(l-s} / A \ " 

Tcr^c'iadl)]^^ i—j (WX)^-^' . (11) 

For s = 0, we get ac(l) = and the exponential behavior 
■m{T,e = 0,x,w) ^ {4wxTy^^^e~*^ . (12) 

If one takes a; = before s — is taken, one gets a 
singular behavior, m{T,e = 0,x = 0,w) — 1/2 for r > 
and to(t, e = 0, x = 0,w) =0 for t — 0. 

To obtain the exponent 6, we take t — 0. In the limit 
e,x,w ^ 0, we get for s ^ 1/2, 



TO(r = 0,e,x,w) = - 



1/3 



(13) 



For < s < 1/2, we get 
to(t = 0, e, X, w) 

-C-^ [ae(l)]~TT7 (ti;x)-Tfj (-S-)T+f , 

The crossover bias in Eq.(14) reads as 

1 + 3 

e 3 / A \ 2(2s-l) 3s 

— = c"[q;c(1)]^(T^ — (wa;) 2(1-2=) 

At s = 0, we have 

TO(r = 0, e, X, w) 



e 

'2A' 



(e < Ecr)] 

(e > 

(14) 

, (15) 
(16) 



m(r,e = 0,x,u;) = 1/ ^/2ac{l)T^''^ ■ 



(9) 



In the equations above, c, c', and c" are constants which 
are independent of x, w, and A. 

These expressions are consistent with the numerical so- 
lution of the mean-field spin-boson model and subsequent 
Aft-scaling analysis (the exponent for x fitted in Ref. [2ll 
deviates due to numerical errors.). A new result here is 
that w also becomes a scaling variable. This implies that 
the logarithmic discretization can no longer keep the uni- 
versal properties intact in the regime ^ s < 1/2 as is 
usually assumed in the NRG studies. 

The above results clearly show that in the regime ^ 
s < 1/2, boson state truncation indeed overtakes the 
critical fixed point of the mean-field spin-boson model. 
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changing it from Gaussian to interacting. For any finite 
Ni,, we get the exponents /3 = (f — s)/(2s) and S — 
(f + s)/(l — s) as long as r or e is sufficiently small. A 
remarkable observation is that these expressions agree 
well with Pnrg and 6mrg, the BNRG exponents for the 
full spin-boson model with finite Nb (See Fig.f4(a) and 
(b))<2i Hence we have the following relations, 

PNRG = {1-S)/{2S) , 

5nrg = {1 + s)/{1~s) . (17) 

In contrast, the exponents Pmf = 1/2 and 5mf — 3 
only appear when t and e are larger than their respective 
crossover scales. 

It is noted that for a finite N^, the result of the mean- 
field Hamiltonian Eq.(5) depends on the parametrization 
used for the bath degrees of freedom. Expressions (9)- 
(16) hold only for H!^Y in Eq.(6), which is obtained from 
a specific parametrization, i.e., the logarithmic discretiza- 
tion. For other parametrization schemes, different expo- 
nents will obtain. At A^f, = oo, bosons become canonical 
and one gets the exact classical exponents irrespective of 
the parametrization of the bath. 

The above discussions are for the mean-field star- 
Hamiltonian obtained from logarithmic discretization. 
The mean-field chain-Hamiltonian cannot be solved an- 
alytically at finite Nt,. Using BNRG, we managed to 
solve the mean field equations by iteration. We found 
that the converged solution has no qualitative difference 
from that of the mean-field star-Hamiltonian. That is, 
the same nonclassical (classical) exponents are obtained 
in the low (high) energy regime. 

B. NRG results for s = 0.3 

In this section, we present the BNRG data for 
m(T, e, A, X, w), carry out scaling analysis, and extract 
the exponents P and 5 in the limit Ni, = oo. The scal- 
ing analysis is in parallel with that in Ref. [2l]. We will 
demonstrate that in the regime < s < 1/2, BNRG data 
fulfills Eqs.(lO) and (14), except that ac(l) there should 
be replaced with the corresponding BNRG values. 

For simplicity, we focus on s = 0.3, a generic value 
in the regime < s < 1/2. Extensive BNRG calcula- 
tions are done with various parameters A^, A, and Ms- 
Due to computational limitations, we use Nf, — 8 ^ 100, 
A = 2 ~ 10, and = 60 - 120. We define t = 
a-ac{A,A,Nb,Ms). ac{A,A,Nb,Ms) is the critical dis- 
sipation strength for a fixed set of parameters. We found 
that ac{A, A, Nb, Ms) has almost no dependence on Nb 
(less than 10~^ percent between Nb = 13 and A'f, = 50), 
similar to the mean-field case where ac strictly does not 
depend on Nb- For A > 2 and A^f, < 100, ac converges 
very fast as Alg increases. For A and A dependence, we 
find 

r Al-^ (A^O); 
ac{A,A,Nb,Ms)x < (18) 
[ ae(l) + c(A-l), (A^l). 



E 



1 

S=0.3 




1 1 




5.212) 






II 






£-4.2 
-^-4.4 




1 




-2.0 -1.6 -1.2 -0.8 



I 1 1 1 : i : — 1_ 

-6 -5 -4 -3 -2 -1 

Ig X 



FIG. 1: lgm(e = 0) vs. Igr at s = 0.3 for different Ab val- 
ues. Throughout this paper, Ig stands for logarithm base 
10. From bottom to top, A^ = 8, 13, 20, 32, 50, 70, 100, re- 
spectively. Dots with solid guiding lines are BNRG data. 
The dashed line is the power law fit for A;, = 100 in small 
r regime, which gives Pnrg ~ 1.169 being consistent with 
(1 — s)/ (2s) = 1.167. The guiding dash-dotted line marks the 
/? = 1/2 behavior in the large A;, limit. Other parameters are 
A = 0.1, A = 10.0, and Ms = 80. Inset: lgm(e = 0) as a 
function of 1/-/Vi, at Igr = —5.212. The guiding dashed line 
with slope —1/2 shows the asymptotic power law in the large 
Nb limit. 

Here c « 0.003 for Nb = 8, Ms = 80 is quite smaU. This 
explains the very good agreement between the BNRG 
curve ac{s) and that from other methodsJ^^— 

1. m(r, e = 0, A, a;, w) 

In Fig.l, we show the Nb dependence of the function 
m{T), with other parameters fixed. We observe that in 
the small t limit, m(T) c>c r'^™"'' and the slope is inde- 
pendent of Nb- To an accuracy of 0.1%, the extracted 
exponent Pnrg agrees with (1 — s)/(2s). We checked 
other s values in the regime < s < 1/2 and the fitted 
Pnrg fulfill Eq.(17) very well, as shown in Fig. 14(a). As 
discussed in Appendix, the expressions in Eq.(17) are the 
results of a specific parametrization for the bath degrees 
of freedom, i.e., the logarithmic discretization. In the 
regime < s < 1/2, different /? and 5 will be produced if 
different parameterization is used. Thus the coincidence 
with Eq.(17) hints that /3 and 5 produced by BNRG are 
actually artificial ones induced by the boson state trun- 
cation and their values dominated by the logarithmic dis- 
cretization. Further evidence of the interplay between Nb 
and A is given below. 

For a fixed Igr, m(r) cx x^^/"^ in the large A^ limit, 
as shown in the inset of Fig.l. Therefore, in the small t 
limit, we have a double power form 

m(T,e = 0, A,a;,w) cx T^"«'^a;-^/^ (19) 
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FIG. 2: Data collapse of lgm(e = 0) vs. Igr curves shown 
in Fig.l. The dashed line has a slope (1 — s)/(2s) = 1.167. 
Inset: the ratios of the shift magnitudes S\gx/S\gT (squares) 
and 6lgm/5\gT (dots) as functions of Nt- The dashed lines 
are a^/ar = (1 — 2s)/s = 1.333 (upper) and Om/ar = 1/2 
(lower), respectively. 
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FIG. 3: g — Igm—am/uAgT vs. z = Igx — a^/arlgT for various 
Nb^s at s = 0.3, A = 0.1 (squares with solid guiding lines). 
Here am/ar ~ 1/2 and a^/ar = (1 — 2s)/s are used. From 
bottom to top, A'^i, = 50, 70, 100, respectively. The dashed 
lines with slope s/(2(l — 2s)) and —1/2 mark the large and 
small r limit, respectively. As A^t increases, the regime with 
almost zero slope expands. 



As N}) increases, the curve shifts along certain directions, 
signaling the scaling behavior. The upper part of the 
curve has an approximate slope Pmf = 1/2. Its range 
is enlarged as iVf, increase. These features resemble what 
was found in the mean-field spin-boson modeli^i It is then 
expected that the crossover Tcr between the two power 
law regimes: the lower one with Pnrg and the upper one 
with I3mf, moves toward zero as x tends to zero. Thus 
at x = the classical exponent ISmf will be recovered in 
the whole r regime. 

Following Ref. [U, we assume that to(t, x) is a general- 
ized homogeneous function (GHF), i.e., m(rA''^,xA°=") — 
\°-"^m(T,x) for any positive A. Letting A = r^/"^, we get 
the scaling form 

fn{f,x) = ^f + -g{x-^f). (20) 

Here, t = logj^gt {t = m, r, x). Or and Ux are the scaling 
powers for r and x, respectively. g{z) is a universal func- 
tion. Using information from Fig.l, i.e., the double power 
form Eq.(19) in small r limit and the non-sigularity of 
m(r, x) in the limit x — > 0, it is easy to obtain the fol- 
lowing ansatz for g{z), 

{const., (Zsat < z < ^cr); 
(21) 
9z, (z>Zcr). 

Here z — x — ax/cirf. In the saturation regime where 
z < Zsat, g oc {am/ax)iz - x). From Eqs.(17),(19), and 
(20), one extracts the exponents 6 = —1/2 and 

am/cLr = 13m F = 1/2 

ax/ar = (1 - 2s)/s . (22) 



One way to verify the scaling ansatz is to show the data 
collapse. The GHF assumption implies that 

m{f + a^X, X + a^X) = m{f, x) + Qm^- (23) 

Therefore, in the log-log diagram, the group of curves 
TO(r, x) should collapse when f , x, and rh are shifted by 
Sf, 5x, and (5m, respectively. The ratios between any 
two of them give the corresponding exponents, 5fh/5f = 
O'm/O'T and Sx/Sf = ax jar- In Fig. 2, a perfect data col- 
lapse is obtained from the data in Fig.l. The ratios 
of the shifts are plotted as functions of iV;, in the in- 
set. Compared with Eq.(22), the agreement is poorer 
for smaller iVf,, probably due to nonuniversal corrections. 
It improves continually as TVfc increases. This forms a 
consistent confirmation of the GHF assumption and the 
results Eqs.(17),(19)-(22). 

The universal function g{z) is plotted in Fig. 3 for N\, > 
50. The downturn on the left part of the curve comes 
from the saturation of m in the large r regime. As x gets 
smaller, the intermediate regime with zero slope extends, 
forming a pronounced plateau as described by Eq.(21). 
In the r7i(r) curve, this corresponds to the extension of 
the regime with /3 — 1/2 a.s Nj, increases. 

The two-section behavior of g{z) in Eq.(21) agrees 
with Eq.(lG). Putting Eq.(21) into Eq.(20), one gets 
m{T, x) cx jj-i/S-r/^NRG for r < Tcr and to(t, x) oc t'^""" 
for r > Tcr. The crossover Tcr is determined by z^r — z^r 
and z — X — ax/arf. We get oc x"^/"^ = x*''^~^^', 
same as in Eq.(ll). 

Guided by the mean-field results Eq.(lO), we also carry 
out scaling analysis for m{T, w) with respect to w = A — 1 
(for s = 0.3 and fixed N),). In Fig. 4(a), similar scaling be- 
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FIG. 4: A scaling behavior of the order parameter at s = 
0.3, A = 0.1. (a) lgm(e = 0) vs. Igr for various A's (Nt = 
8, Ms = 80). Symbols with guiding lines are BNRG data 
and dashed lines are fitted lines. Inset: lgm(e — 0) at Igr — 
—5.212 as functions of Igw for Ni, = 8, Ms = 80 (squares) and 
Nt = 20, Ms = 120 (dots), respectively. The slopes of the 
fitting solid lines are —0.46. Other parameters are the same 
as in Fig.l. (b) Universal function of g — \gm — am / aAgT v.s. 
z — Igw — aw/aAgT for different A's. Here am/or = 1/2 and 
a^/ar = (1 — 2s) /s are used. The dashed lines with slope 
s/(2(l — 2s)) and —1/2 mark the large and small r limit, 
respectively. As w decreases, the regime with almost zero 
slope expands. 



havior as in m(T, x) is observed. In the inset, Igm is plot- 
ted as a function of \gw for a fixed Igr, giving a power law 
behavior with exponent —0.46, consistent with the —1/2 
in Eq.(lO) within numerical errors. By assuming that 
TO(r, w) is a GHF with scaUng exponent {am, o-t, o,w), us- 
ing the exponents Eq. (22) , and by applying the data col- 
lapse procedure (not shown), we obtain a^/ar = 1/2 and 
Q-io/q-t = (1 ^ 2s)/ s — axjobT- In Fig-4(b), the universal 
function concerning TO(r, w) is plotted using the above 
exponents. Similar to Fig. 3, the universal function has 
the form of Eq.(21), with x ^ w. Fig. 4 shows that the 
scaling variable w plays a similar role as x. Therefore, the 
crossover scale Tct has an additional factor ty^/f^^^s) p^j, 
either w = or a; = 0, the TO(r, a;, vj) curve will have the 
classical exponent /3mf in the regime r <C T^at- Here, 
again the NRG data agree with the mean-field expres- 
sions (10) and (11). The fact that w becomes a scaling 
variable means the failure of the logarithmic discretiza- 
tion in NRG: A may alter the universal quantities. This 
is solely due to the boson state truncation. As learned 
from the mean-field study, the parametrization scheme 
for the bath is relevant for the critical behavior once the 
bosons are no longer canonical. 




-5 I • ' ' ' ' ' ' ' ' ' • 1 

■6 -5 -4 -3 -2 -1 

IgA 



FIG. 5: (a) lgm(e = 0) vs. Igr for different A's at s = 
0.3. From bottom to top, A = 10"\ 10"^, 10"^, 10""*, 10"^, 
respectively. Symbols are BNRG data and the dashed lines 
are power law fit from small r data, which agree with Pmrg = 
(1 — s)/(2s). Inset: Igm at Igr = —16 as a function of A 
(symbols). The fitted line (solid line) has a slope —0.821 
which agrees well with —(1 — s)^/(2s) = —0.817. (b)lgac as a 
function of IgA. Symbols are BNRG data and the solid line 
is the fitted line with slope 0.703 which agrees with 1 — s. 
Calculated with Nt = 8, = 80, and A = 4.0. 

As for A, we show the BNRG results in Fig. 5. The 
power law fits in Fig. 5 (a) and its inset disclose a double 
power behavior in the small r limit: m{T, A) oc A'^t^"'^'^ . 
Here 77 is an exponent to be determined below. In 
Fig.5(b), we show that BNRG produces ac(l) oc A^~^, 
a result that has been established by BNRG and pertur- 
bative renormalization group study.— To understand the 
exponent 77, we resort to Eq.(lO). If we collect the ex- 
ponent of A in Eq.(lO), and are careful enough to use 
the BNRG result adl) = c(A/a;c)^"'' (instead of the 
mean-field one ac(l) = As/(2wc)) in Eq.(lO), we get 
77 = -(1 - s)V(2s) = -0.817 for s = 0.3. The exponent 
fitted from BNRG data is —0.821, as shown in the inset 
of Fig. 5 (a). The excellent agreement supports that the 
A-dependence in BNRG results can also be summarized 
by Eq.(lO). We checked other s values in < s < 1 and 
confirmed this conclusion. 

Combining the above results, we draw the conclusion 
that in the regime < s < 1/2, BNRG resutls for 
m{T, e = 0, A, w, x) are well described by Eqs.(lO), except 
that the A-dependence of ac(l) in the equation should be 
replaced by the BNRG form Eq.(18). Correspondingly, 
the crossover point Tcr reads as 

cx Ai-"(wa;)"/(i-2s)^ (24) 
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FIG. 6: lgm(r = 0) vs. Ige for different iV^'s at s = 0.3, A = 
0.1. From bottom to top, Nt = 8,13,20,32,50,100, respec- 
tively. Other parameters are A = 0.1, A = 10, Ms = 80. 
Symbols with solid guiding lines are BNRG data. The dashed 
line is power law fit for Ni, = 100 in small e limit, which gives 
slope 1/5nrg = 0.533 being consistent with (1 — s)/(l + s). 
The guiding dash-dottod hno marks the asymptotic behav- 
ior 1/(5 — 1/3 in the large Nt hmit. Inset: lgm(T = 0) at 
e = 10^^ as functions of Iga; for A = 10, Ms = 80 (scjuares) 
and A = 4, Ms = 120 (dots), respectively. The dashed lines 
are guiding lines with slope — s/(l + s) = —0.231 to mark the 
power law behavior in large Nb limit. 



As a scaling variable, A is different from w or x. As 
seen in Fig. 5(a), when A decreases, both the Pnrg- and 
/3MF-exponent regimes in m(T) curve shift to the left 
along the horizontal direction, but the ^MF-exponent 
regime does not expand. This is because another 
crossover scale T,at, above which m(r) ^ 1/2, also de- 
creases as A^~*. Indeed, using to(t, A) oc [T/ac(I)]^/^ oc 
(r/Ai-^)V2 in Tcr <^ T <^ Tsat and m(r,A) ~ 1/2 in 

T > Tsat, one gets Tsat OC A^~*, Same as Tcr- 



2. m(r = 0, £, A, a, w) 

In this part, we fix r = and study the order parame- 
ter m as a function of e. This is related to the exponent 
5. 

In Fig. 6, m(e, x) versus e at t = is plotted for var- 
ious X values. Similar to m(T,x), in the small-e limit 
all curves for different A^b's have the same power law 
m{e,x) cx 6^/*""° with x-dependent coefficients. For 
s — 0.3, the fitted power 0.533 agrees well with the mean- 
field expression (I — s)/(l-|-s) = 0.538. For other s values 
in < s < 1/2, the comparison is shown in Fig. 14(b). 
The inset of Fig. 6 supports m(e,a;) oc xP for fixed e and 
large A^b, with the exponent Q' to be determined below. 
Therefore, a double power form in the small e limit is 
obtained, TO(e,a;) oc x^ fr(i-'')/(i+s) i^^^ other limit 
where e is much larger, the upper dashed line in Fig. 6 



marks out a finite regime where TO(e, x) oc e^/'^mf — gi/3 
There is a crossover ecr oc which 
separates the lower double-power regime from the upper 
classical regime, goes to zero as x — > 0. 

To obtain 0', we carry out the scaling analysis based on 
GHF and do the data collapse in Fig. 7. We assume that 
ra{e,x) is a GHF with scaling powers (a^, 0^,03;) and 
get m{e,x) = e°-^l'^^h{xlt°"'l°-^). The universal function 
h{z) has a two-regime form similar to Eq.(21), i.e., 



const., {Zsat <C ^ < Zcr); 

h{z) oc { (25) 

d'z, [z » Zcr)- 

Here z = x— (oa;/oe)e. Using m(e, 0) oc e^/^ and am/dx = 
s/[2(l — 2s)] confirmed before, we obtain 

^-s/{i+s)^\/5mrg^ [e<^ecr); 
m(e, a;) oc ^ (26) 



.1/3 



(e > tcr)- 



Here ecr oc x^^/'^". This means 6' = — s/(l + s) and 



Qm _ 1 

tte 3' 

_ 2(1 - 2s) 
a. 3s 



(27) 



In the inset of Fig. 6, the BNRG result for 6' (symbols) 

compares favorably with the exponent — s/(l+s) (dashed 
lines). In the inset of Fig. 7, the exponents am/ae and 
Ux/ae extracted from ratios of the shifts 6m, de, and 
Sx agree well with Eq.(27) in the large Nb limit, con- 
firming Eq.(26). The similarity to the mean- field equa- 
tion Eq.(14j is obvious. Using this 9' we also obtain 
Ccr oc a;^^/[^'-^~^*)l. It agrees with the mean-field result in 
Eq.(15). 

The A scaling analysis is carried out in Fig. 8 for 
m(T = 0,e,w) for fixed A, A^;, and Ms. In Fig.8(a), 
the small e regime of the curves fulfills the power law 
m(e, w) oc ei/''«RG! ^nd the larger e regime m(e, w) oc 
^ However, for a fixed e, the shift of the curve with 
A is very weak for A between 2.5 and 10. In Fig. 8(b), the 
order parameter m(e = 10~ w) does not fit a power law 
of w as nicely as m(T = to,w) does in Fig. 4. The data in 
the small A regime are Mg dependent. Here the crudest 
estimation from larger A regime gives the power —0.07. 
Assuming a GHF form for m(e, w) (at fixed r = and 
a; > 0) and using previously obtained relations among 
Um jCLr^O'e and ttw , we get 



m(e, w;) oc e 



1/Snrg,„-s/(1+s) 



(28) 



same as in the mean-field result Eq.(14). The fitted 
power —0.07 is quite difi'erent from — s/ (l + s) = —0.231. 

We believe that this discrepancy is due to the fact that 
for 2 ^ A ^ 10, w is not small enough to enter the scaling 
regime for observing the correct power law. This is in 
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-4 -2 

Ig e (shifted) 

FIG. 7: Data collapse of lgm(r = 0) vs. Ige curves shown in 
Fig.6. The dashed straight line has a slope (1 — s)/(l + s) = 
0.538. Inset: the ratio of the shift magnitudes 5\gx/5lge 
(squares), and Slgrn/Slge (dots) as functions of Ni,. The 
dashed straight lines are a^/a^ = 2(1 — 2s)/(3s) = 0.889 
(upper) and a-m/ae = 1/3 (lower), respectively. 




for different A's at 
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= 0.3. 
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FIG. 9: lgm(T = 0) vs. Ige 

From bottom to top, A = 10"^ 10"^, 10"^, 10"", and 10 
respectively. Symbols are BNRG data and dashed lines are 
power law fit, giving an average slope 0.539 which agrees with 
(1 - s)/(l + s)= 0.538. Inset: lgm(T = 0) at Ige = -11.292 
as a function of IgA (symbols). The fitting (solid line) gives 
a slope —0.534, being consistent with e/A scaling. Other 
parameters are Ni, = 8, and Ms = 80. 




FIG. 8: (a)lgm(r = 0) vs. Ige for different A's at s = 0.3. 
(b)lgm(r = 0) vs. Igw at e = 10"'^ for various (iV^, M<,)'s. The 
fitted slope is about —0.07 in the shown regime, different from 
— s/(l + s) = —0.231. This is probably due to the fact that w 
is not small enough. Inset: the mean-field results Igm vs. Igw 
for T = 0, e = lO"'^ (cycles) and r = 10~^,e = (squares), 
respectively. Other parameters are A = 0.1, Nb = 8. 



contrast to Fig. 4 where a nice w-scaling persists up to 

A = 10. We have not fully understood this contrast yet, 
but only mention that we observed similar differences in 



the mean-field calculations. There, the shift with w is 

much more pronounced in m{T = To,e = 0,w) than in 
m(r = 0, e = cqj w). The curves Igm — \gw of the mean- 
field Hamiltonian are plotted in the inset of Fig. 8(b), for 
T = (cycles) and e = (squares), respectively. It is 
seen that the expected power law behavior appears only 
when A <C 2. In the regime 2 < A < 10, the lgm(T = 
0,e = eo,w) curve significantly deviates from the correct 
power law behavior, resembling the BNRG results. This 
supports om notion that much smaller A is required to 
observe the m{T = 0,e = eQ,w) oc behavior. 
As a consequence of this reasoning, w should also appear 
in the crossover scale as a factor 

We checked our results for Nb and Ms up to 50 and 
300, respectively, and found that the quality of the data 
is not improved. In the BNRG calculations, after each 
diagonalization, only the lowest 1/Nh fraction of eigen- 
states are kept. Hence, as Ni, increases, one needs larger 
A or larger Mg to compensate the error from discarding 
states. As a result, it is very difficult to produce reliable 
data in the large-A'^;, and small-A regimes. For a fixed 
A^b, it is known that to approach a smaller A regime, one 
needs to use larger Mg. However, we did not find the 
expected power law behavior of w up to A = 2 using 
ATfe = 8 and = 300. 

In Fig. 9, we study the A scahng in m{e, A) at r = 0. 
The figure resembles that of lgm(r, A) — IgA. The double 
power form in the small e regime is found to fulfill the 
e/A scaling, i.e., 



m(e, A) (X 



(29) 



with 1/Snrg = (1 ^ *)/(l + s). As A decreases, 
both I/Smf- and l/cJAr/jG-exponent regimes move toward 
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smaller e along the horizontal direction. Simple analy- 
sis shows that a crossover scale €sat oc A separates the 
saturation regime m ^1/2 from the l/i5MF-exponent 
regime. Ccr must have the same factor A because the 
classical exponent regime is not enlarged as A deceases. 
Summarizing these results, we get Csat > ^cr and 



Esat OC A . 



(30) 



These results concerning the scaling behavior with A is 
also consistent with the mean-field expression Eq.(14)- 
(15), provided that we use the BNRG result ac(l) oc 
A^~* in the equations. 



3. Summary for < s < 1/2 

We summarize the above analysis. Our main conclu- 
sion is that in the regime < s < 1/2, due to the bo- 
son state truncation, the order parameter m produced by 
BNRG is a scaling function of variables t, e, A, a; = l/A^h, 
and m; = A — 1. Interestingly, this function agrees well 
with the mean-field equations for finite N^, [Eqs.(lO) and 
(11) and (14) and (15)], except that the A-dependence 
of q;c(1) should be replaced by the corresponding BNRG 
one, i.e., ac(l) oc A^"**. In the small-r or -e limits, the 
boson-state truncation A^f, introduces artificial exponents 
B and 5. which arc different from the correct classical val- 
ues. We summarize the BNRG results inO<s< 1/2 as 
the following: 



m(T, e = 0, A,x, w) 

A"Vr5, (-Tcr- < r «; Tsat); 

I 1/2, (r > Tsat), 



= < 



with the two crossover scales Tsat > Tcr given by 

Tcr oc A^'"" {WX)T^ , 

For m(r = 0, e, A, x, w) we have 



m(r = 0,e,x,w) 

-cA~T+j(y;a;)-TTJ(-^)Tqrf ^ (c <C Ec 



(31) 



(32) 



c'(i)V3, 



A 

I -1/2, 



(ccr < e < esot); 
(e > £sat)- 

(33) 
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FIG. 10: lgm(e = 0) vs. Igr at s = 0.7. (a) A = 10.0 
and Nt = 8, 13, 20, 32, 50,rcspoctivoly; (b) A^t = 8 and A = 
3.0, 4.0, 6.0, 10.0, respectively. Symbols are BNRG data and 
the dashed lines are power law fits. Other parameters are 
A = 0.1, and Ms = 80. 



with the two crossover scales Csat > Ccr given by 



Ccr OC {WX) 2(1-23) A , 
Csat OC A . 



(34) 



For the role of A in the order parameter to, we 
observed t/A^~* scaling in the function to(t, e = 
0,A,w,x), while the function to(t = Q,e,A,w,x) has 
e/A scaling. These features are independent of the bo- 
son state truncation, and hence should hold also in the 
regime l/2<s<l. As we will see in the next section, 
this is indeed the case. 



C. BNRG results for s = 0.7 

In this section, to compare with the < s < 1/2 case, 
we study s = 0.7, a generic value in the regime 1/2 < s < 
1. In this regime, the mean-field theory predicts classical 
exponents (3=1/2 and 5 = 3 for any Nb- The boson 
state truncation does not influence the Gaussian critical 
fixed point in the mean-field Hamiltonian. 

For the full spin-boson model, BNRG predicts an in- 
teracting critical fixed point and nonclassical exponents 
/? and 5. In Figs. 10(a) and 10(b), lgTO(r, e = 0) — Igr 
curves are plotted for various Nb^s (A = 10) and A's 
(Nb = 8), respectively. Other parameters arc fixed. It 
is clearly seen that the situation is dramatically different 
from s = 0.3: there is no Nb or A scaling. Therefore, the 
boson state truncation and A do not influence the correct 
extraction of (3. 
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FIG. 11: lgm(e = 0) vs. Igr at s = 0.7 for different A's. 
From bottom to top, A = 10~\ 10~^ 10"^, 10"*, 10"^, re- 
spectively. Symbols are BNRG data and the dashed lines are 
power law fit, giving an average slope 0.296. Inset on top 
left: IgOc vs. IgA data gives the fitted slope 0.307, consistent 
with 1 — s. Inset on bottom right: Igm at Igr — —6.949 as a 
function of IgA (symbols). The fitted slope is —0.096, close to 
—f]{l — s) = —0.089, as expected from t/A}'" scaling. Other 
parameters are A = 4.0, Nt, — 8, Ms = 80. 
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FIG. 13: lgm(r = 0) vs. Ige at s = 0.7 for different A's. 
Symbols are BNRG data and dashed lines are power law fit. 
From bottom to top, A = 10"\ 10"^, 10"^, 10"*, 10"^, re- 
spectively. Inset: lgm(r = 0) at Ige = —10.423 as a function 
of IgA. The solid line is a power law fit in small A limit and 
the slope -0.179 is close to -(1 - s)/(l + s) = -0.177, as 
expected from e/A scaling. Other parameters are A = 4.0, 
iVft = 8, and Ms = 80. 



scaling, as observed in the s = 0.3 case, it is then easy to 
obtain 
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FIG. 12: lgm(r = 0) vs. Ige at s = 0.7 for (a) different iV^'s 
and A = 10.0, and (b) different A's and A^;, = 8. Symbols are 
BNRG data and the dashed lines are power law fit, giving the 
average slopes 0.175 and 0.173 for (a) and (b), respectively. 
Other parameters are A = 0.1 and Ms — 80. 



m(r,e = 0,A)cx A-(i-^)'3t'^, {r-^T^at). (35) 

For s = 0.7, the fitted value of rj is —0.096, which agrees 
with —(1 — s)f3 ~ —0.089 quite well. We also checked 
other s values in the regime 1/2 < s < 1 and found good 
agreement, confirming the validity of r/A^"* scaling in 
this case. The relation ac(l) oc A^"*^ is demonstrated 
in the top left inset of Fig. 11. Tsat oc A^~* is the 
crossover scale separating the power-law regime from the 
saturation regime m 1/2. 

Similar analysis is carried out for m(T = 0, e, A, x, w). 
In Figs. 12(a) and 12(b), we plot this function for s = 0.7 
and A 0.1, for various A^b's (A = 10) and A's (iVf, = 8). 
No Nf, or A scaling is observed, similar to the curves 
m(T, e = 0, A, X, w) shown in Fig. 10. 

The A scaling carried out in Fig. 13 shows that the ex- 
ponent 6 fulfills the expression (l-|-s)/(l — s), even in the 
1/2 < s < 1 regime. This is consistent with the previous 
findings of the hyperscaling relation S = (1 + x)/{l — x) 
and X = Si^ In the inset of Fig. 13, Igm — IgA is shown for 
fixed e, disclosing a power law behavior consistent with 
m(T = 0,e,A) oc A^'-^-'''^/'-^+'\ This implies the e/A 
scaling in the regime 1/2 < s < 1. Therefore, we have 



In Fig.ll, we study the A scaling at s — 0.7 for 
m{T,e = 0). The main plot and the bottom right in- 
set show that m{T,e — 0,A) oc A^t^. In the regime 
1/2 < s < 1, exact expression for /3 as a function of s is 
not known. If we assume that m(T,e,A) has (r/A^-*) 



m{T = 0,€, A, w,w) oc y-^y^" (e<esat)- (36) 

Here Esat esc A is the crossover scale separating the power 
law regime from the saturation regime m 1/2. 
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D. 



s = and s = 1/2 



The boundary cases s — and s — 1/2 need some 
discussions. For s — 0, the mean-field solution at fi- 
nite Nb gives [Eq.(12)] m{T, e = 0) cx exp(— A/4u;cr) and 
the crossover scale Tcr becomes independent of wx. This 
means that is large and decreases with wx at sublead- 
ing order. We carried out the BNRG calculation for s = 
with 10-^ sC A < 10^2 and 10 A^b «C 50. Using the ex- 
ponential fit, we obtain Uc — within an error less than 
10^^. For a fixed Ni,, m{T, e = 0) shows a nice exponen- 
tial behavior m oc exp(— c/t) with c « (0.27 ± 0.02) A. 
This is consistent with the mean-field expression (12). 
For the prefactor, we observe the x"^/^ behavior using 
8 < A^b ^ 50, and find the —1/2 exponent of w using 
A < 2.5 and Ms > 200. However, we also find a per- 
fect t/A scaling at s = 0, being different from Eq. (12). 
By combining these results, we expect the crossover scale 
Tcr cx: —A/\n{wx). As for m(T = 0, e), BNRG gives the 
exact result m{T — 0,e) = — e/(2A) since at r = a = 
the spin is decoupled from the bath. 

For s = 1/2, contrary to s = 0, the powers of wx in 
the expressions for Tcr and ecr diverge, leading to the 
vanishing of the latter for wx < 1. Therefore, the clas- 
sical regimes in m(T, e = 0) and m(T = 0, e) expand to 
the zero r or e limit, even for finite N^s and A > 1. 
Indeed, using Ni, — 8 and A = 4, BNRG gives accurate 
/?MF and 6mf at s = 1/2, as shown in Fig. 14(a) and (b), 
respectively. For s = 1/2, corrections to power laws in 
m(T, e = 0) and m(T = 0, e) may arise in the subleading 
order in which x and w may play some roles. However, 
such effects are difficult to observe in the BNRG data 
because of numerical errors and fitting errors. 



III. SUMMARY AND DISCUSSION 

In Fig. 14, we summarize the exponents /?, 6, 7 and x 
for the spin-boson model in the whole regime < s < 1. 
The squares in the figure show the naive BNRG results 
obtained using finite Nj,. The solid lines for (3 and 6 
in the regime < s < 1/2, as shown in Fig. 14(a) and 
(b), are ^ = 2s/(l - s) and S = {1 + s)/(l - s), the 
exact solution of the mean-field Hamiltonian with finite 
Ni,. They agree well with the BNRG data (squares). 
The dashed lines are the correct values extracted from 
the scaling analysis. They deviate significantly from the 
BNRG results in ^ s < 1/2. In the regime 1/2 ^ s < 1, 
the boson state truncation error and mass fiow error still 
exist but they do not influence the exponents. 

We also studied the magnetic susceptibility x using 
BNRG in the regime < s ^ 1/2 (squares in Fig.l4(c)). 
It is found that xijiT — 0) curve is independent of w 
and X and 7 = 1 holds at high precision. For the A de- 
pendence, BNRG calculations for various s in the regime 
^ s < 1 confirm the scaling form x c< A~^(t/A^~'*)~^, 
being consistent with m(r, e. A) = /(t/A^~^, e/A). In 
the mean-field theory, the factors containing wx cancel 
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FIG. 14: Critical exponents of the spin-boson model as func- 
tions of s: (a) /3, (b) 5, (c) 7, and (d) x. The squares are 
the exponents fitted from BNRG data with finite A'';,. The 
solid lines are Pnrg = (1 - s)/(2s), &nrg = (1 + s)/(l — s), 
'JNRG = 1, and xnrg = s, respectively. The dashed lines are 
the correct exponents free from boson state truncation error 
and the mass flow error. 



each other and one obtains x ^ ^^^i'''/'^c{^))~^^ inde- 
pendent of wx. It covers the BNRG result when ac{l) 
is replaced by A^~^. In the regime 1/2 < s < 1, the 
BNRG result for 7 follows a nonclassical curve [squares 
and the dashed line in Fig. 14(c)], the exact expression of 
which is not known yet. It is related to /3 and S through 
hyperscaling relations. 

The temperature dependence x(T) at t = is gov- 
erned by the exponent x. We did not find any observable 
shift in x(t = 0, T) with the changing of Nb or A in the 
range 8 < A^fc < 20 and 2 < A < 10. We always obtain 
a; = s in the regime < s < 1/2. This indicates that 
the mass-fiow problem is deeply rooted in the NRG algo- 
rithm and can not be solved by the scaling method that 
we discuss in this paper. According to the analysis of the 
mass-flow problem in NRG, x — 1/2 should hold in 
the regime ^ s ^ 1/2 and a; = sinl/2<s<l. This 
conclusion together with the BNRG result are summa- 
rized in Fig. 14(d) for completeness. 

To conclude, we confirm the statements in Ref. [l^ 
that for the spin-boson model, the critical exponents /3 
and S are classical in the regime ^ s ^ 1/2, while 
are nonclassical inl/2<s< 1. BNRG with finite 
Nb produces artificial interacting critical fixed point in 
the regime ^ s < 1/2, and the extracted /3 and S 
are incorrect. In the regime 1/2 ^ s < 1, boson-state 
truncation does not play a role and BNRG results are 
correct. Scaling analysis with respect to x ~ 1/Nb and 
w = A ~ I can be used as a supplement to ordinary 
BNRG to overcome the boson-state truncation error. It 
is also observed that the order parameter has the scaling 
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form m(r, e, A) = /{t/A^ ^^e/A) in the whole regime 
^ s < 1. 

Two features of the BNRG results can be understood 
in terms of the mean-field Hamiltonian. First, as shown 
in Eqs. (31) and (33), the boson-state truncation error 
disappears in the limit A = 1. Second, the truncation Ni, 
influences only the exponents but not the critical cou- 
pling ac- In the mean-field Hamiltonian [Eqs. (6)-(8)], 
the nth boson mode has a displacement ^n/S,n and its 
coupling strength to the spin is proportional to 7„ . Both 
7„ and ^„ decay exponentially with increasing n, but the 
ratio 7n/^n diverges as A^^"*'"/^. The low-energy modes 
(with large n) are more displaced and need more states 
to describe. Therefore, truncation of the state infiuences 
mostly the low-energy modes and not the high-energy 
modes. In the limit A = 1, however, there is no diver- 
gence in the displacement and all the modes can be de- 
scribed equally well by a finite N},, hence the truncation 
error disappears. Since the critical exponents are closely 
related to the energy distribution and the displacement 
behavior of the low-energy boson modes, they are sus- 
ceptible to the truncation. While the critical coupling 
is dominated by the integral properties of the spectrum 
in which the high-energy modes have more weights, it is 
thus only very weakly influenced by the truncation. 

Our results support the validity of quantum-to- 
classical mapping for the spin-boson model in the regime 
^ s < 1/2. Note that our picture of m(T) for the trun- 
cated spin-boson model is consistent with the results of 
exact diagonalization study in Ref. [l^. There, although 
a diff'erent basis set is used, in principle, the boson-state 
truncation error also exists and the nonclassical expo- 
nents should be present in the low-energy limit. However, 
the fitting of data relatively far away from the critical 
point (note the linear scale in Fig. 3 of Ref. 15) coinci- 
dentally misses the nonclassical regime and produces the 
correct classical exponent 13m f- 

For other quantum impurity models such as the 
Bose-Fermi-Kondo model (BFKM) model~22 and its 
anisotropic version, the Ising-BFKM,^^ similar quantum- 
to-classical mapping arguments supports Gaussian fixed 
points in the regime ^ s < 1/2. Correspondingly, one 
expects /3 = l/2,(5 = 3 and x — 1/2. Neither hyper- 
scaling relations nor uj/T scaling should hold. However, 
variouls studies using BNRG^^, QMC simulatioi*S^"— , e- 
expansion^^'^^ and the large-N limit analysis2^i21i2^ all 
point to the failure of the quantum-to-classical map- 
ping in these systems: they obtained interacting fixed 
point and nonclassical exponents. Especially, for the 
anisotropic BFKM with a sub-Ohmic boson bath, the 
exponents are found to be the same as l3pfB.G and Snfi.g 
of the spin-boson model. 

In light of the present study, it is possible that the 
same boson state truncation error may exist in the BNRG 
study of the Ising BFKM in Ref. [13 and lead to incorrect 
/? and S in the regime < s < 1/2. Note that the boson 
state truncation problem cannot be remedied by going 
to smaller energy scaleS i ---'"- Since the fermion bath in 



the Ising BFKM can be integrated out to produce an 
additional boson bath with Ohmic spectrum, the critical 
behavior should be dominated by the sub-Ohmic one, 
leading to the expectation that it belongs to the same 
universality class as the spin-boson model. In the quan- 
tum Monte Carlo simulations, it was observed that the 
truncation of correlations in the imaginary time axis is 
the key to produce the interacting critical point»2i Other- 
wise Gaussian behavior will obtain. These may be signa- 
tures that the quantum-to-classical mapping holds also 
for the Ising BFKM in < s < 1/2. It is straightforward 
to examine this statement using BNRG supplemented 
with the Nf, scaling. 

As far as the isotropic BFKM is concerned, the sit- 
uation seems different. Here, the symmetry is different 
from the spin-boson model and the Berry phase effect 
is claimed to be nontriviali^ Much work has been done 
for this model, supporting the failure of the quantum-to- 
classical mapping in the sub-Ohmic regime.— Due 
to the very subtle nature of this problem, however, ex- 
act numerical studies are still desirable as a further 
confirmation. Both the isotropic and the anisotropic 
BFKM are at the core of studying the Kondo lattice 
with Ruderman-Kittel-Kasuya-Yoshida (RKKY) inter- 
actions using the extended dynamical mean-field the- 
Qj.y ^25, 32^34, 35^37- _39 Therefore, any definite conclusions 
concerning these impurity models will have important 
impact on the understanding of the competition between 
Kondo screening and the antiferromagnetic state in the 
heavy-fermions metals. 

In summary, we carried out systematic BNRG studies 
for the spin-boson model in the sub-Ohmic regime, sup- 
plemented with the scaling analysis for the boson state 
truncation x — l/Ni, and the logarithmic discretization 
parameter w = A — 1. For < s < 1/2, the func- 
tion m{T,€ = 0,A,x,w) [m(r = 0, e, A,x,u')] is shown 
to bear a multiple power form in the small-r (-e) limit. 
Classical exponent /3 = 1/2 (5 = 3) is identified in the 
regime r 3> (e 3> fcr), agreeing with the conclusion 
from the quantum-to-classical mapping. The crossover 
scale Tcr i^cr) gocs to zero in the small x or w lim- 
its in a power law. This presents a scenario of how 
the boson-state truncation error disappears in the limit 
Nh — > cxD. The observation that x and w always ap- 
pear as a product (wx) indicates that, in the regime 
^ s < 1/2, the boson-state truncation invalidates the 
logarithmic discretization scheme, which is the basis of 
NRG. Independent of the issue of boson-state truncation, 
we also find that the scaling form for the order param- 
eter m(T, e. A) — /(r/A-'^^'^, e/A) in the whole regime 
«C s < 1. 
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Appendix A: Critical Exponents of the Mean-Field 
Spin-Boson Model 

In this appendix, we present the calculation of the crit- 
ical exponents /3, S, and 7 for the mean-field spin-boson 
Hamiltonian. 

The Hamiltonian of the spin-boson model reads as 



H 



sb 



(Al) 



The mean-field Hamiltonian is (neglecting a constant) 



1 + Jl]A.(a. + al) 



For a finite iVb, the critical behavior depends severely on 
the parameterization scheme for {uji,\i\ as well as on 
the discretization scheme. In order to compare the expo- 
nents with the NRG ones, we carry out the logarithmic 
discretization as in NRG and get 



(ct ) °° 



where the superscript LD denotes the logarithmic dis- 
cretization. 7„ and ^„ are obtained from standard proce- 
dure*'" and given in Eq.(7) and (8) in the main text. In 
Eq.(A3), the spin and boson degrees of freedom decouple 
and we get Hj^'j = Hspm + Hboson, where 



Hspin — ^l^^x + + d)crz 



and 



00 ^ oc 

Hboson = E ^nOl^^n + ^(^2) ^ ("« 



n=0 



n=0 



Here 



ri=0 ^ 



(A4) 



(A5) 



(A6) 



For Nb ^ 00, Hboson can be solved exactly. For finite 
Nb, it is no longer exactly solvable. In Ref. [2l|, H^^^ is 
solved numerically and the critical behavior is studied. 
Here, we start from the single mode Hamiltonian 



(A7) 



In the small 6 limit, the ground state is not influenced 
by the truncation and we have (a^ -I- 0)0 — —26 and 
(a^a)G — 9^. In the large 9 limit, the number of bosons 
in the ground state is confined by the boson state trun- 
cation Nb- Therefore in this limit (a^a)G = Nb. For the 
displaced harmonic osc illator, to leading order of Nb we 
have {a + a^)G = —2a/ {a'^ajc = — 2\/]Vb. The crossover 
scale is determined by equating the two situations. In 
summary, we have 



[a' +a)G = 



^29, 
'2^b, 



(A8) 



» ^Nb) 



We confirmed this result numerically. Here {■■■)g denotes 
the ground state average. {al^+an)G has different form in 
the regimes n ^ no and n ^ no. Here no separates the 
freely biased mode (small n) from the saturated biased 
mode (large n). We can then solve Hboson approximately 
and obtain one of the mean-field equations from Eq.(A6), 



n=0 



in 



n— no + l 



1ni 



(A9) 



with no determined by equating ^/Nb with the effective 
bias, 



7no 



1, 



= \ Nb 



(AlO) 



Introducing x = 1/Nb and m = (cr^) /2 and summing over 
n in Eq.(A9) produces 

d = ai(A)aTO + [a2(A) + &i(A)] aT^xT^m^ . (All) 
The parameters ai(A), a2(A), and foi(A) read 

4(2-1-5) [l-A-(i+^)]' 



ai(A) 



(1 + 5)2 [i_ A-«] [1- A-(2+'^)] 



c^c, (A12) 



a2(A) = 2^(l-fs)-^(2-f s)^ 



1 - A"(^+'^) 



1 - A-(i+-^) 



A'' - 1 



and 



61(A) = -a2(A) — 



A* - 1 



(A13) 



(A14) 



'A(i+s)/2 _ 1' 

The other mean-field equation from solving Hspin is 

2 

(A15) 



1 



A2 



m = — 
2 



+ v/A2 + (e + d)2 



A2 



e + d+ ^A^ + {e + dy 



2 • 



Near the critical point, it reduces to 



e + d 1 /e + d 
" 2A ^ 4 I A 



O 



e + d 



(A16) 
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Putting Eq.A(ll) into the above equation and keeping 
only the leading term of each type, one gets 



m : 



e ai{A)am 
2A 2A 
a2(A) + bi(A) 
2A 



+ 



4A3 



[ai(A)am]' 



l+£ 



(A17) 



From this equation the critical coupling strength is ob- 
tained as ac(A) = — 2A/ai(A). It is remarkable that it is 
independent of the boson state truncation Nb. Introduc- 
ing T = a — ac{A.) to measure the distance to the critical 
point, we get the self-consistent equation as 



e a2(A) + 61(A) 



[ac(A)] xi-»mi- 



2t . % 

H TTrm - 4m = 0. 

ac(A) 



(A18) 



To investigate the role of A, we define w = A — 1 and 
expand adA), a2(A) and 61(A) to leading order of w. 
Finally we get the self-consistent equation in terms of r, 
e. A, w, and x as (keeping ac{A) in the definition of r 
unchanged) 

2t , q e r / s _2_ i±£ 

— —m - 4m'^ - -r - c— ac(l) (wx) 1-= m i-« =0. 
Q!c(l) A A 

(A19) 

Here ac(l) = A.')y(2co'c) is the critical a value for A = 
1. c = 2(2-'')/(i-'') [l/s - 2/(1 + s)] is a constant. It 
is noted here that the anomalous term involving the 
{wxyl^^~^^ comes solely from the factor a2(A) -|- 61(A), 
instead of from expansions of etc (A) in Eq.(A18). Solv- 
ing this equation in various limits produces the analytical 
expressions Eqs.(9)-(16). 

Eq.(A19) reduces to standard mean- field equation and 
gives /3 = l/2,(5 = 3, when the rrt' term dominates over 
the anomalous term (^wxYl'^^~''^'rrS^^^^^'^^^^'' . This is the 
case when m is small in the regime l/2<s<l, or when 
wx <C m in the regime < s < 1/2. The anomalous term 
will dominate over the w? term and give /? = (1 — s)/ (2,s) 
and (5 = (1 -|- s)/(l — s), when m — >■ in the regime 
< s < 1/2. Considering different regimes of r and e, 
we get the expressions (9) and (10) and (13) and (14). 
The crossover scales Tcr [Eq. (11)] and tcr [Eq.(15)] are 
obtained by equating the expressions in the classical and 
nonclassical regimes. 

At s = 0, Eq. A(10) is replaced with 



no 



1 

inX 



In 



(l + A-i)2 1 



1 - A-i 8aa;m2 



(A20) 



The summation in Eq.A(9) produces 

1-A-i , r A(l 



d = —Samujc 
—SarriLOr 



(l-^A-i)lnA 
A-1/2 + A-1 

1-FA-i • 



In 



+ A-i)2 



1 



A ^)xam'^ 

(A21) 



This gives c^c = 0. Expanding at A = 1 to leading order 
and combining Eq.(A16), we get the order parameter for 
s = 0, 



1 _1 A_ 

m(r, e = 0, A, a;, w) = c-(wxt) , (A22) 



and 



m(T = 0, e. A, x,w) = — 



2 A' 



(A23) 



Here, c is a constant independent of A, x, and w. Note 
that if one takes the a; ^ limit first and then takes the 
limit s 0, Eq.(A22) becomes 



1/2, (r > 0); 
m{T,e = 0,A,x = 0,w) = { (A24) 

0, (r = 0). 



The magnetic susceptibility x can be obtained by tak- 
ing derivative on both sides of Eq. (A19) with respect to 
e. This leads to the expression 



-Ax = 

C(l + S) U>c r /,M^^ / ^ s 2s O 2t 

-i; '--^ ac 1) {wx) — m— + 12m^ 

1-s A «c(lj 



1 -1 



(A25) 



Here, c = 2(2-^)/(i-'*) [l/s - 2/(1 -h s)]. Analyzing this 
equation in different t regimes, we obtain x(t, A) = 
— ac(l)/(4AT) in the regime 1/2 ^ s < 1, and 



=(1) 



x(t,A) 



4At 



: OcO) 



4s At 



(t » Tcr)-, 
(r < Tcr). 



(A26) 



in the regime ^ s < 1/2. This gives the exponent 

7 = 1, independent of Nh and w. 

We also studied the critical behavior of m(T, e, A, x) 
using other parametrization and discretization schemes 
for the bath spectrum. It is found that for finite Nb, the 
critical exponents l3 and 6 are strongly dependent on the 
scheme. For Nb = 00, they all reduce to the mean-field 
values. 
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